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Recent experiments report the current (I) versus 
voltage (V) characteristics of a tunnel junction 
consisting of a metallic tip placed at a distance d 
from a planar electrode, d varying over six orders of 
magnitude, from few nanometres to few millimetres. 
In the 'electric-field-assisted' (or 'field emission') 
regime, as opposed to the direct tunnelling regime 
used in conventional scanning tunnelling microscopy, 
all I-V curves are found to collapse onto one single 
graph when d is suitably rescaled, suggesting that 
the current I = l(V,d) is in reality a generalized 
homogeneous function of one single variable, i.e. 
I = l(V-d- k ), where A. being some characteristic 
exponent and T{x) being a scaling function. In this 
paper, we provide a comprehensive explanation — 
based on analytical arguments, numerical simulations 
and further experimental results — for the scaling 
behaviour that we show to emerge for a variety of 
tip-plane geometries and thus seems to be a general 
feature of electric-field-assisted tunnelling. 



1. Introduction 

A sharp tip approached vertically to a conducting surface 
at subnanometre distances and biased with a small 
voltage with respect to the surface builds a junction 
across which electrons can be transferred from the tip 
apex to the nearest surface atom (or vice versa) by 
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Figure 1. Schematic of the tunnel junction. The emitter (shaded in grey) possesses rotational symmetry with respect to the 
z-axis and carries thex, y A z-coordinate system at its apex. The azimuthal angle 0 used throughout the paper is defined as 
the angle formed by the x-axis and the projection onto the Ay-plane of the coordinate vector of an observation point in space 
(black dot). In a typical l-V experiment, a positive voltage V is applied to the planar counterelectrode residing at a distance d 
from the emitter and the resulting current / flowing from the planar counterelectrode to the emitter (meaning that the electrons 
flow from the emitter to the planar counterelectrode) is recorded. 

direct quantum-mechanical tunnelling. Such a junction is used, e.g., in scanning tunnelling 
microscopy (STM) for imaging the surface topography with the spectacular atomic spatial 
resolution that was awarded the Nobel Prize in 1986 [1-3]. When the distance d between 
tip and collector is increased, one enters the electric-field-assisted tunnelling regime [1,3,4], 
where the current is dominated by electrons emitted from the (typically sharp) tip into the 
vacuum region residing between the tip apex and the target through a classically forbidden 
zone enveloping the tip apex (figure 1). Such a regime is, for instance, the one underlying 
the topografiner technology [5,6] — an imaging technique which was the precursor of STM 
but was abandoned, probably because of the enormous success of STM. Electric-field-assisted 
quantum tunnelling is also widely used in recent and less recent developments in micro- and 
nano-electronics [7,8]. Recent experiments [9] in the regime of electric-field-assisted tunnelling 
suggest a remarkable scaling invariance of the current flow with respect to changes in the 
tip-to-collector distance d by several orders of magnitude (from a few nanometres to a 
few millimetres). This scaling invariance was detected by observing the collapsing of the 
family of l-V curves, measured at various distances d, onto one single curve when the 
voltage was suitably rescaled with a scaling factor R depending on d through a power law 
R ~ d~ x , i.e. I = X(V • d~ x ), X(x) being a scaling function. Such a scaling invariance — well known, 
e.g., in the field of critical phenomena [10] — is certainly not realized in the direct tunnelling 
regime [1], and it is also not usual in solid-state electronics, so that its observation is yet 
somewhat surprising and unexplained. It is the scope of this paper to provide a comprehensive 
explanation of the experimentally reported [9] scaling invariance by using analytical and 
numerical arguments, and by introducing further experimental results that establish the scaling 
behaviour as a systematic property of junctions in the electric-field-assisted tunnelling regime. 
The paper is organized as follows. 

In §2, we consider the electrostatic problem of a sharp metallic tip, forming one side of the 
junction, which is approached vertically at a distance d by a planar counterelectrode, forming 
the other side of the junction. We formally solve the associated Laplace equation and find 
analytical expressions for the electrostatic potential <P{x,y,z) by considering tip and plane as 
equipotential boundaries, with the tip being at ground and the plane being at a positive potential 
+V. Within our treatment, we consider only highly symmetric, realistic tip shapes. The wording 
'highly symmetric realistic 7 describes here the fact that the geometries considered in §2 are (i) close 
to the shape one expects for 'real' tips (as revealed by a systematic tip imaging via light and 
electron microscopy [9]) and (ii) sufficiently symmetric so that the electrostatic problem can be 
solved to a large extent analytically. On the one hand, this choice helps the reader to follow the 



main arguments without relying on numerical results only, and on the other hand it is well- 
known that unequivocally revealing power laws in experiments and numerical computations 
is a very difficult task, but it is quite straightforward within an analytical framework. The 
main result of §2 is the conjecture that in the vicinity of the tip end the 'almost' singularity 
associated to the tip apex introduces a power-law dependence of the potential on d and V of 
the form: 

<P(x,y,z;d,V)~V - d~ x • <P(x,y,z), (1.1) 

where A is a characteristic exponent that depends on the geometrical details of the tip and on 
the range of distances d (see the detailed discussion in §2). <t>{x,y,z) is some function containing 
only the coordinates (x, y, z), (x, y, z) indicating points residing very close to the tip apex so that 
z<^d. The meaning of the symbol '~' in equation (1.1) and other equations in this paper will be 
discussed in detail in §2. Note that, for example, for a 'planar' tip (in which case we should no 
longer speak of 'sharp tip', of course) X assumes the trivial value of 1. 

In §3, we use the one-dimensional Jeffreys-Wentzel-Kramers-Brillouin approximation of 
quantum tunnelling [11] to compute the tunnelling current density / from the tip into the planar 
electrode in the presence of a potential that behaves as given in equation (1.1). The main result 
of this section is to show that, for an emitter subject to the electrostatic potential worked out for 
different geometries in §2, the scaling behaviour of the potential implies that the current density 
/ becomes a function of one single variable 

]{ y fd) = J{ y. d -^ f (1.2) 

where X is the same exponent appearing in equation (1.1). We obtain an analytical expression 
for the scaling function J{x) both for X = 1 [12-15] and for These exact scaling results led 

us to identify a fundamental length scale A<p — the De Broglie wavelength associated with the 
maximum height of the tunnelling barrier (cp) — and an effective barrier width A that determine 
the leading behaviour of the current density 

J(A (p ,A)~e- A / A *. (1.3) 

This equation is certainly true for the few, highly symmetric models of electric-field-assisted 
tunnelling discussed in this paper, but it might hold approximately in general and for both direct 
and electric-field-assisted tunnelling. 

In §4, we present experimental data and numerical results from simulations that corroborate 
the scaling hypothesis and establish its systematicness. Moreover, we report on the experimental 
observation of small deviations from the scaling hypothesis and discuss their possible origin. 
Finally, we attempt to combine experiments, numerical results and scaling hypothesis into a 
proposal for a functional dependence J(x) that explains most of the data. 

The appendices A and B present the mathematical details of the calculations leading to the 
results summarized in the main text. The scope of these appendices is to allow the verification 
of our computations which, reporting about the relatively recent idea of scaling in the tunnelling 
regime (Cabrera et al. [9]), might evoke some (salutary) scepticism. 

2. Electrostatics of the junction in the electric-field-assisted tunnelling regime 

Within the purpose of this paper, we consider a conducting tip as an infinitely long object with 
'small' cross section and ending with a more or less sharp apex. The tip is kept at 0 potential and it 
is placed vertically at a distance d from a conducting plane held at potential + V (figure 1). The aim 
of this section is to evidence the scaling behaviour of the potential <t> (x, y, z) on the parameters V 
and d, which are typically imposed experimentally [9]. If Q denotes the region of space excluding 
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Figure 2. (a) Emitter with conical profile. The coordinates of a point (black dot) in the y — z-plane are described by the 
intersection of a line of constant spherical coordinate # and a circle of radius r centred at the origin. The full angle of aperture of 
the cone co Q is indicated in the figure, [b) Emitter with cuspidal profile. The coordinates of a point (black dot) in they — z-plane 
are described by the intersection of a circle of constant bispherical coordinate (circle centred at the focal point / = a) and a 
circle of constant bispherical coordinate a (a circle with centre along the axis/ = 0). Tangent lines to the cusp at the tip apex 
define an angle co Ql which can be identified as the full angle of aperture of the cusp. 



the tip and the plane, then the electrostatic problem defining the electrostatic potential 0 (x, y, z) is 
a well-defined Dirichlet problem and reads 

V 2 <2> = 0 in^2 

0=0 on the surface of the tip 

0 = +@ d on the plane 

and \<P{x,y,z)\ < +0^ V(x,y,z)eG, 

where the last equation follows from the maximum principle of harmonic functions. The solution 
of equation (2.1) is unique and can be computed, at least numerically, in the entire space Q. 
However, for the purposes of §§3 and 4, where analytical expressions for the tunnelling current 
density will be derived, only the behaviour of the potential in the very vicinity of the tip apex and 
along the tip axis is required [1,11]. In this section, we therefore focus only on the behaviour of 
0 (x = 0, y = 0,z) for small z. In the following, we summarize the main results. The details of the 
calculations are presented in appendix A. 

Remark. The boundary condition on the tip assigns a uniform value of the electrostatic 
potential to the tip surface. We use the symbol V (figure 1) to denote the experimental 
voltage recorded during the measurements of the experimental I-V characteristics, instead. We 
acknowledge that V contains an electrostatic contribution (which can be considered to be the 0j 
defined here), but also a contribution due to work function difference between tip and planar 
electrode, not considered in these sections. Moreover, real emitters might have a non-uniform 
spatial distribution of work functions, which is also neglected in this section. We point out that 
the difference between 0^ and V is typically 1 V or less, so that the distinction between them can 
usually be neglected. However, for small d and for small V, this difference might be a source 
of small deviations from practical scaling (experimental deviations are indeed observed and 
reported later). 



(a) Conical and cuspidal tips 

These geometries describe the shape of real tips on the micrometre scale [9], but contain a true 
singularity at the tip apex (figure 2). For both distant (d » 0) and near (d small) planes, we have 
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Figure 3. (a) Small sphere of radius a terminating a conical profile with aperture angle co 0 . [b) Hyperboloid of revolution. A 
point (black dot) in the/ — y-plane is the intersection of a line of constant prolate spheroidal coordinate v (an ellipse with focal 
points at ±a along the axis/) and a line of constant spheroidal coordinate v (a hyperbola with focus at / = tfor/ = — a). 
The asymptotes of the hyperboloidal profile define the full angle of aperture of the tip ap- 



proved in ref. [9] (see also appendix A for more details) that the potential in the vicinity of the tip 
apex has a leading term of the form 

<P(z,d)~<P d -(?f\ (2.2) 

where the exponent X\ is defined as the smallest index X for which the Legendre function Px(x) has 
a zero at x = cos(tt — &>o/2). For small angles of aperture, the exponent k\ is given approximately 
by [16,17] k^coo) = [21n(2/^ 0 )] _1 , whereas for coq = 7i , corresponding to a planar emitter, we have 
\\ = 1 [16,17]. 

Remark. Quoting from ref. [18, pp. 39-40], almost Verbatim' but 'mutatis mutandis', we 
point out that the relation <P{z,d) ~ &d • (z/d) kl does not imply the relation <P{z,d) = A • (z/d) Xl . 
In general, we find (see ref. [9] and appendix A) that there is an infinite number of additional 
correction terms of the type (z/d) Xk , with — and this is the crucial point — > X\. In the limit 
(z/d) -> 0 — which is the asymptotic case considered for discussing 'scaling' — these additional 
terms become infinitesimally small with respect to the leading one (z/d) Xl , so that we can write 

0(z,d) 

lim \ — ~r = const., (2.3) 

z/d^O 0 d • (z/d)^ 

which is equivalent to writing &(z,d) ~ <S>& • (z/d) kl . Equation (2.3) establishes the significance of 
the symbol '~' and our use of it in relation to the concept of 'scaling'. 



(b) Sphere-on-the-cone and hyperboloid of revolution 

The sphere-on-the-cone model is particularly suitable for mimicking a rounded tip with overall 
conical shape (figure 3). In fact, the sphere-on-the-cone model terminates the cone with a small 
sphere of radius a. Hyperboloids of revolution are suitable to mimic rounded tips with asymptotic 
conical shape. They have, in fact, two asymptotes that can be used to define a full angle of aperture 
coq. The two asymptotes meet at a point in front of the apex that is the intersection between the 
z^axis and the so-called confocal plane. The focal length a is the distance between this point and 
the focal point of the hyperboloid, which is located within the tip on the tip axis. In contrast with 
conical and cuspidal tips, the scaling properties of the leading potential term close to the tip apex 
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Figure 4. Paraboloidal emitter. A point (black dot) in the z — y-plane is the intersection of a line of constant paraboloidal 
coordinate u (a downward oriented parabola) and a line of constant paraboloidal coordinate v (an upward oriented parabola). 
The tip apex is at a distance /? 0 /2 from the origin of the coordinate system, where R Q is the radius of curvature of the tip. 
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for the sphere-on-the-cone and the hyperboloid of revolution depend on whether the plane is 
'distant' (d > a) or 'near' (d<&a). In particular, we obtain 



for d ^ a 



/a\^i z 

and <P{z, d) ~ &d • - • - for d > a, 

\d/ a 

where k\ is the same exponent found for conical and cuspidal tips. 

(c) Paraboloid of revolution 

A paraboloid of revolution is characterized by a radius of curvature Ro, but it does not have a 
characteristic angle of aperture and therefore misses one essential characteristic of the tips used, 
for example in ref. [9] (figure 4). However, it provides an interesting limiting case 



(2.4) 



and 



®{z,d)> 
(P{z,d)' 



ln(2d/jRo) 



z 

Ro 



for d <$C Ro 
for d^>R 0 . 



(2.5) 



In fact, by exploiting the identity lim^oW*^ — 1)/a0 — mx / we rnight consider the logarithmic 
dependence on d in equation (2.5) as a special case of equation (1.1) when X -> 0. 

(d) General tip geometry 

Equations (2.2), (2.4) and (2.5) provide analytical expressions for the potential along the tip axis, 
ready to be used for computing the tunnelling current (§3). However, the scaling behaviour with 
&d and d observed along the tip axis can be extended by continuity to a small neighbourhood of 
the tip apex as well, yielding equation (1.1). Moreover, the scaling hypothesis, equation (1.1), has 
been verified explicitly for a restricted number of highly symmetric geometries, but we propose 
that it might have general validity. 
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. The tunnelling current density in the presence of a non-trivial exponent k] 

One of the most remarkable results of §2 is that, even if the planar electrode is a 'distant' one, the 
boundary condition on the plane determines the electrostatic potential in the vicinity of the tip 
apex, where electric-field-assisted (or field emission) quantum tunnelling occurs. Having shown 
some scaling properties of the electrostatic potential with d, now we would like to address the 
question whether these scaling properties affect the field emission process at all. In this section, 
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Figure 5. Sketch of the one-dimensional tunnelling barrier profile along the tip axis, showing the potential energy originating 
from the electrostatic potential for the case of a linear dependence on z, for simplicity. This sketch is used to define some 
quantities used in the bulk of the section. Z| and z 2 are the points at which the potential energy — e<P(z) crosses the Fermi 
level, with |z 2 — z^ \ being the width of the potential barrier. <p is the work function of the tip. 



we compute the tunnelling current density for a very simplified model of electric-field-assisted 
quantum tunnelling, which only involves the presence of the pure electrostatic potential within 
the tunnelling barrier (figure 5). We also assume, for the sake of simplicity, the purely one- 
dimensional model used, for example, in ref. [1], which foresees free electrons within the tip, 
and is therefore strictly applicable only to a large flat planar emitter [19]. For very small radii 
of curvature — in particular for the conical and cuspidal models — quantum confinement can, 
for example, occur and the electron energies might become quantized [20]. We also neglect the 
image potential correction to the electrostatic potential energy, which is known to substantially 
lower the barrier height and to modify by orders of magnitude the current density [1,11]. Thus, 
owing to these assumptions, the calculations performed in this section are far from being realistic. 
Nevertheless, we would like to point out that the scaling properties derived in this section 
on the base of the pure electrostatic potential are actually obeyed by the experimental data 
presented in ref. [9] and in this paper (with the limitations discussed in %3a). Given this agreement 
with experiments, we suspect (but cannot prove it!) that all the elements neglected by our 
simple — perhaps trivial, model — are not really modifying the scaling properties obtained from 
the pure electrostatic potential. Of course, we are convinced that this last sentence will raise some 
controversial discussion, but we feel that it might justify the publication of the results presented 
in this section. 

Within our simple, not-so-realistic-model, the current density along the tip axis can be 
written as [1] 



\Jz\=e 



dv z • v z • p(v z ) • D{v z ), 



(3.1) 



where — e is the charge of the electron, v z is the z-component of the velocity of the free 
electrons within the metal and dv z p(v z ) is the number of electrons per unit volume with velocity 
between v z and v z + dv z . Finally, D(v z ) is the transmission coefficient of tunnelling describing 
the probability that an electron with velocity v z overcomes the potential barrier. By inserting the 
suitable (free electron like) expression for p(v z ) [1,11] in equation (3.1) and by using the standard 
Jeffreys- Wentzel-Kramers-Brillouin approximation of quantum tunnelling [1,11], one obtains the 
following result: 



\h\=h(E { ).e- G ^ 



G(E f )= 2. 



l87T 2 m 



\j-e0{z)-E i dz 



(3.2) 



where Ef is the Fermi energy of the tip, — e0 (z) is the potential energy of the electron within the 
tunnelling barrier, m is the mass of the electron and G(Ef ) is the so-called Gamov exponent. The 
integration limits Z\ and z^_ are the zeros of the quantity — e&(z) — Ef (figure 5). 

We now summarize the main scaling results obtained for the current density by assuming a 
potential that satisfies the scaling hypothesis derived in the previous section. The details of the 
derivation are described in appendix B. 



These two geometries for the tip provide a special case. First, the apex represents a geometrical 
singularity where it is difficult to imagine that a 'current density' is well defined, as the 
area of the apex itself is zero. Second, the electrostatic potential is non-analytic within the 
tunnelling barrier: &(z) ~ 0& • (z/d) Xl . The electric field at the apex, e.g., diverges to infinity, so 
that standard formulae of field emission [1,11,14,15] are useless. We are therefore in the presence 
of a mathematically difficult problem: 'a zero area' of emission and an 'infinite' electric field. It 
is therefore most remarkable that application of equation (3.2) to the non-analytic electrostatic 
potential (note that the integrals in equation (3.2) can be computed exactly in terms of known 
mathematical functions (see equations (B3) and (B4))) produces a finite current density ]{0^,d) 
which has a remarkable scaling property associated with the scaling behaviour of the electrostatic 
potential: although there are two independent experimental variables 0^ and d, the current 
density is a function of one single variable 0& • d~ Xl , namely 



where a\,a2 are dimensionless numbers, containing, for example, natural constants and J coryQ {x) 
is the scaling function. The • d~ Xl scaling agrees with the experimentally observed scaling 
behaviour [9]. The scaling function will be discussed in §4. 

Remark. It is beyond the scope of this paper to attempt to derive a more realistic expression 
for the emission current from a conical or cuspidal emitter. The significance of equation (3.3) 
is that, even in a situation of extreme singularity, the simple model underlying equation (3.2) 
provides the base for the data collapsing observed experimentally. Note that the calculation 
has produced a fundamental scale for the energy, cp, and a fundamental scale for the length 
A(p = y/h 2 /2m(p. Typical values [9] for these parameters are cp ~ 4.5 eV and A v & 0.6 nm. A v is the 
De Broglie wavelength corresponding to the maximum barrier height <p. Finally, we note that the 
case corresponding to a planar electron emitter in front of a planar counterelectrode [13-15] is 
also covered by the scaling law when the suitable value k\ = 1 is inserted. 



Realistic tips — such as those described by the sphere-on-the-cone, hyperboloid and parabolic 
models — might have some rounding, characterized by a spatial scale that we called a in §2 and is 
absent in the conical geometry. Accordingly, as shown in §2, the rounding produces for all of them 
a potential with leading linear term within the tunnelling barrier. The scaling of the potential with 
d depends on whether the planar counterelectrode providing one boundary condition is distant 
(0(z) ~ (a/d) kl • 0^ • z/a) or near (0(z) ~ 0& • z/d). In both cases, the integrals in equation (3.2) 
can also be computed in terms of elementary functions [11,13] and the resulting scaling 



(a) Non-analytic potential: conical and cuspidal tip 



and 




(3.3) 



(b) Rounded tips: analytical potential 



laws write 



i/ distant (^d^)i = j' distant * e0d ( a 



(3.4) 



and J dis ^' near (x) = af stant * near • * 2 • c^ - "* 

Note that for near planes, the scale a cancels out from the problem and only one fundamental scale 
remains, namely A v . By contrast, for distant planes, the problem has two fundamental scales: a — 
originating from the geometry of the tip — and A<p — originating within the quantum mechanics of 
the tunnelling process. The dimensionless constants a\ and depend on the geometry, on natural 
constants and, in the case of distant planes, on a/A^. 

(c) Proposal for a unified approach and relationship with direct tunnelling 

If we introduce the barrier length A = |Z2 — z\\, the tunnelling current densities for the non- 
analytic potential (equations (B3) and (B4)) and for the linear potential (equation (B7)) can be 
unified into one single expression in the variables <p (the barrier height) and A: 

7(A^) = hf •^•r^H (3.5) 

where b\ and b^ are some numbers that take into account the exact shape of the potential within 
the tunnelling barrier. Note that the exponential (leading) term in equation (3.5) appears also 
in models of direct quantum tunnelling, the process underlying STM (see eqn. (25) in ref. [1]), 
where A is given by the tip-to-surface distance. Therefore, it seems that equation (3.5) — at 
least its exponential factor — is a general motive in tunnelling phenomena. What distinguishes 
direct tunnelling from electric-field-assisted tunnelling — the process described in this paper — is 
the functional dependence of / on <&<x,d. In the case of electric-assisted tunnelling, the scaling 
hypothesis imposes the special functional dependence of A on <Pd,d, i.e. A = A(@d • d~ x ). The 
current density behaves accordingly and the functional dependence is reduced to one single 
variable. In direct tunnelling, this reduction to one single variable is not realized (see eqn (25) 
in ref. [1] and ref. [3]). We also mention here that the spin of the electrons is a further variable that 
appears in direct (STM) tunnelling phenomena [21]. It would be interesting to study whether the 
spin also plays a role in electric-field-assisted tunnelling and its scaling properties. 



4. Experimental and numerical evidence of the scaling hypothesis 

(a) Experimental evidence 

(i) /-I/ curves 

The first observation of the scaling behaviour is contained in the experimental data reported in 
ref. [9]. We reproduce part of these data in figure 6fl(i,ii), which also contains further experimental 
data (fr(i,ii) an d c(i,ii)). Most remarkably, these further data, taken with different tips, strongly 
support the collapsing behaviour reported in ref. [9]. We conclude that the scaling behaviour is a 
general feature of electric-field-assisted tunnelling. Let us discuss now some details of the data in 
figure 6. The top panel of the figure shows a set of current (I) versus voltage (V) curves taken for 
three different tips (a,b,c) at different distances from a doped Si(lll)-single crystal surface used 
as counterelectrode. Experimental details about the measurements are reported in ref. [9]. We 
recall that our tips are fabricated starting with a tungsten wire with a few millimetres length and 
250 \im diameter. The last few hundreds of micrometres close to one end of the wire are etched 
electrochemically to assume a cuspidal profile which, in the final few micrometres towards the 
apex, resembles very much a cone with a full angle of aperture between 6° and 12°. Electron 
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Figure 6. {a-c)(\) Family of l-V curves in the range of d = 3-300 nm. [a-c)(\\) The curves on the top are made to collapse 
onto a reference curve (chosen arbitrarily as the one corresponding to d = 300 nm) by multiplying the voltage with a number 
R[d), plotted in the insets on the bottom. The continuous curves through R[d) in the insets are power laws a c/~\ with 
X « 0.2 ± 0.05. The alternative horizontal scale in the bottom figure gives the electric field F at the apex, as derived in ref. [9]. 



microscope imaging of the tip reveals a rounding of the tip towards its apex. The rounding varies 
between 5 and 30 nm, depending on the details of the tip preparation in ultra-high vacuum [9]. 
The counterelectrode is typically a W(110) or a Si(lll) single-crystal surface: it turned out that 
low-noise l-V characteristics are favoured by the use of typically very flat Si(lll)-surfaces. The 
family of curves in the top of the figure collapses onto one single reference curve (figure 6a-c{\i)) 
when the voltage is multiplied by a ^-dependent factor R{d), which is well described by a power 
law of the type ^d~ x , with X ~ 0.2. The analysis of the data presented in figure 6 demonstrates 
that J is a function of V • d~ x , i.e. I = X(V • d~ x ). This is the essential point of the scaling hypothesis 
proposed in §3. Furthermore, the experimentally observed values for the exponent X fall in the 
range of values expected for X\ by equations (2.2) and (2.4). We recall that the scaling hypothesis 
in §3 refers to the current density while experiments measure the total tunnelling current. Yet, the 
scaling hypothesis is realized for the current as well. This points to the fact that the details of 
the area on the apex where the electrons originate from (uniformity of emission, variable size of 
the area of emission with d and V) are not relevant for the scaling behaviour. In other words, the 
tunnelling process is dominated by the exponential function over any multiplicative, probably 
non-exponential prefactor. 

(ii) V-d curves 

The scaling equation I = X(V • d~ x ) can be inverted to yield V = X^ _1 ^(J) • d x , where X^ -1 ) indicates 
the inverse scaling function. Accordingly, if we plot the experimental data of figure 6 in a 
V-d diagram, a family of V-d curves appears (figure 7^(i-iii)). When V is multiplied by the 
scaling factor R(I) = 1/X^~ l \l) (see insets in figure 7fr(i-iii)), the family of V-d curves collapses 
onto one single reference curve, behaving as d x , with the exponent being the same entering 
the insets of figure 6. We find that the factor R(I) closely follows a curve of the type ^J _/x , 
with a very small value for \i [9]. It follows that \/X^~ l \l) ^ J _/x and therefore X(x) ^x 1///x . We 
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Figure 7. tf(i-iii) Family of V-d curves (data taken from figure 6), each curve corresponding to a given current (colour coded 
along the vertical bar), b(i-iii) The family of curves is collapsed onto one single reference curve by multiplying the voltage with 
a number /?(/), plotted in the inset. 



note that, in the range available experimentally, this scaling function describes properly the 
experimental results but is unusual in tunnelling phenomena. On the other side (see ref. [18, 
pp. 40-41] for a complete discussion), a very small critical exponent might be the sign of a peculiar 
singular behaviour, like, for example, a logarithmic singularity. In fact, for small jjl and moderate 
currents, one has J _/x ~ 1 — \i In J, and thus 1/Z( -1 )(J) ^ 1 — \i In I. Inverting the latter, we obtain 
X(x) ~ e -1 /^*. This is the scaling function predicted by equation (3.4) and, in the available range 
of experimental values of I, it appears to describe as properly as the scaling function x 1 ^ the bulk 
of the experimental data. Note, however, that on the basis of the experimental data alone, it is 
impossible to discriminate between the scaling functions x 1///x and e~^^ x . The two functions could 
hypothetically be distinguished in the limit of x going to infinity. However, this limit corresponds 
to very large currents, while a maximal current of only a few microamperes can be driven through 
realistic tips without destroying them. Note also that there are still some details of the reference 
graphs in figures 6 and 7 which are not properly covered by either scaling functions. 



(iii) Deviations from scaling 

In order to enhance these details, we single out two typical I-V curves within a so-called 
Fowler-Nordheim [13] plot of log(J/V 2 ) versus 1/V (figure 8a). Although the scaling function 
of equation (3.4) predicts a strict linearity of the graphs in this kind of plot, we observe a 
systematic downward curvature of the graphs, which is more pronounced at larger distances. The 
observed curvature represents a deviation from the scaling function predicted by equation (3.4), 
while the separation of the two curves on figure 8a for small voltages means a slight deviation 
from collapsing towards smaller currents. In the following, we offer some arguments that might 
account for these observations. 

The scaling function for a non-analytical potential implies naturally the experimentally 
observed downward curvature, but it cannot be applied straightforwardly to real tips, which 
have necessarily a rounded apex. A simple way of avoiding the unphysical non-analyticity by 
simultaneously keeping its welcome nonlinearity is to introduce a finite cut-off length r into the 
cone solution by substituting the power-law (z/d) Xl with its Taylor series at a finite distance r 
from the tip apex. Here, r might be considered a measure for the radius of curvature. This yields 
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Figure 8. (0) Two /-I/ curves selected from figure 6 are plotted in a log 10 (//l/ 2 ) versus 1/1/ plot, showing a small downward 
curvature and a deviation from scaling at low currents, (b) l-V curve measured at d = 175 nm in a junction with W[W) single 
crystal plane as a counterelectrode. The tip used in this experiment is the one described in figure 9. The continuous curve is a fit 
with a function of the type log 10 (//l/ 2 ) = a\ + a 2 • 1/1/ + 03 ■ 1/l/ 2 ,0/ being the fitted parameters. Inset: the l-V data in a 
linear plot. (Online version in colour.) 



where A being a dimensionless factor (see appendix A). Note that a similar quadratic correction 
to the linear term can also be found using the sphere-on-the-cone or the hyperboloidal models 
for a rounded tip, but the mathematics entailed by the physically plausible cut-off approximation 
is more transparent. We also observe that this approximation introduces a linear and a quadratic 
term into the potential within the tunnelling barrier, but the presence of the quadratic term does 
not break the scaling invariance. The simplest way of taking the quadratic term into account in 
the calculation of the tunnelling current density is to use the proposed unifying equation, which 
implies finding the value of z for which the argument of the square root in equation (4.1) vanishes 
(<P(r) is set for convenience to zero by modifying the boundary condition so that the zero of 
the potential is at z = r). Mathematically speaking, this means solving a quadratic equation. In 
the range of V for which the equation has a solution, we approximately find 

A^r ? +r l -^( ? V (42) 



and accordingly 

J(V,d)^exp 



r • (p/kxeVAir/d)^ + r • (1 - A,i)/4(<p /XieVA(r/d) Xl ) 2 
-a, 



(4.3) 



The quadratic part of the potential, which mimics the upward curvature of the potential energy 
entailed by the original power law, increases the tunnelling width above the value obtained 
with the linear term only and introduces an extra term which is responsible for the downward 
curvature observed in the experimental Fowler-Nordheim plots (see figure Sa and ref. [9]). Note 
that equation (4.3) still obeys the scaling hypothesis which is apparently 'robust' also with respect 
to quadratic corrections of the potential. However, the appearance of the non-leading power 
(r/d) 2kl in the expression for the tunnel width suggests that, ultimately, a realistic expression for 
the current density will need the use of other non-leading powers of r/d as well. But as soon as 
we allow non-leading powers of the type (r/d) Xk , the scaling behaviour is broken and deviation 
from collapsing is expected (and indeed observed, see ref. [9] and figure 8a). We conclude that 
a scaling function of the type J{x) = a\ e ~ a 2/x-a 3 /x ^ w j l j c ] 1 generalizes the expression obtained in 
equation (4.3) for a nonlinear potential, might be useful to interpolate experimental data that show 
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Figure 9. (a) Potential profile <P (z) along the tip axis, for a given distance d. The tip used in this simulation is a hyperboloidal 
model of a 'real' tip (a = 1528 nm, co Q = 11°). The planar counterelectrode is moved between d = 20 nm and d = 1646 nm. 
All profiles are made to approximately collapse onto the reference curve (d = 1646 nm) when the potential 0 (z) is multiplied 
by R& (d) (see inset). R& (d) is shown in the top right corner of the figure. Note that both R<p (d) and the scaled profile in the 
inset behave as a power law, with a power-law exponent of about 0.3. This value is larger than the one of about 0.21 expected 
from coq = 11° . We point out that the value of 0.21 is an asymptotic value expected when the planar counterelectrode is much 
farther away than the confocal plane. In the present simulation, d is smaller or in the range of the confocal plane distance so 
that we might not have reached yet the true power law range. In fact, we have performed computations with d = 6000 nm 
(not shown) and the exponent is seen to converge towards the analytical value of 0.21. We note, however, that the collapsing 
of profiles is realized also in this non-asymptotic range, showing that the scaling property itself is a robust one. {b) Scaling of 
0 (z) and @did)- The continuous line is the potential profile obtained as described above, for d = 1646 nm. The full squares are 
0d{d) data points obtained at a given current, the current of 150 pA having being chosen so that the @ d {d) curve almost lies 
onto the <P (z) graph, without need of a rescaling factor. Inset: <P d (d) data points obtained in a junction with Si(111) at selected 
currents, given in the legend. In the inset, data are rescaled so that they fall onto the same power law. The potential profile 
computed for d = 300 nm (continuous curve) can also be rescaled onto the same curve as the @(j[d) data. 



a downward curvature and small deviations from scaling in Fowler-Nordheim plots. An example 
of such a procedure is reported here in figure 8b. A similar proposal was also put forward in 
ref . [22] on the basis of numerical results (see also ref . [23] for a recent discussion on the curvature 
problem). 

(b) Numerical evidence 

The principle of Saint- Venant [24] implies a further interesting scaling symmetry of the electric- 
field-assisted tunnelling junction. In fact, this principle can be used to extend the validity of the 
conical solution, equation (2.2), to real tips, supposing that they can be viewed as a 'cone with a 
rounded apex'. According to the principle of Saint- Venant, if the rounding of the cone singularity 
is local enough — say limited to a scale length q, which does not need to be atomic, describing the 
'radius of curvature' — then the conical solution, equation (2.2), can be used in the range q <^z^d 
as well (the origin of the z-axis being, as usual, the apex of the tip): 

0(z) ~ 0 d • 1 , a « z « d. (4.4) 

This equation is valid for any real tip shape with arbitrary but sufficiently localized rounding, 
provided that asymptotically away from the rounded apex the conical shape is recovered. Note 
that the region of validity of this equation makes it irrelevant for the field emission process, 
which occurs at positive z <^ q, so that the discussion of this paragraph refers only to a property 
of the electrostatic potential. Equation (4.4) provides a further scaling law which can be tested 



numerically by verifying, e.g., the collapsing of the family of <i>(z,rf) -curves taken for different 
d onto one single curve, when <P{z) is multiplied by a suitable R<p{d)^d Xl (figure 9 a). The 
latter is inversely proportional to R(d) introduced previously (figure 6), i.e. R<p{d) ~ 1/R(d). 
For the computation of the potential, we have fitted a hyperboloid of revolution onto the 
electron microscope micrograph of the tip used for the taking of the experimental V-d-curves. 
Subsequently, we have used a numerical routine [25] to compute <P{z) in the presence of a 
plane placed at a well-defined (large) fixed distance in front of the hyperboloidal tip. A further 
consequence of equation (4.4) is that an experiment where d is changed and is adjusted 
so that the tunnelling current (i.e. the potential within the tunnelling barrier) is kept constant 
yields <t>& ~ d k , i.e. the same power-law dependence that one expects for the z-dependence of the 
potential 0 itself. In other words, we expect that all experimental — d graphs (i) can be made 
to collapse onto themselves (figure 7 and inset of figure 9b) and (ii) they can be collapsed onto a 
0(z) profile, provided that one is not too close to z = d or z = 0 (figure 9b). 
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Appendix A. Electrostatics 



We will consider tips ('needles') with 'realistic' but highly symmetric geometries. The most 
important symmetry, common to all tips, is the rotational symmetry with respect to the tip 
axis. One speaks of axially symmetric Laplace problem because the sought for potential will be 
independent on the rotational angle 0 (figure 1) and that part of the Laplace operator which 
contains derivatives with respect to 0 can be neglected from the beginning. A further symmetry 
is imposed on the electrostatic problem: the surface of the tip has a particular geometry along 
which, when suitable, typically curvilinear coordinates x and y are introduced and one of them 
(say, y) assumes a fixed value xjq. One of the boundary conditions reads therefore: &(x,yo) = 0 Vx. 
In view of the need to fulfil this boundary condition, it is suitable to introduce a separation Ansatz 
0(x,y) = X(x) • Y(y) for solving the associated Laplace equation, yielding 



AxX(x) 
X(x) 



AyY(y) 
' Y(Y) ' 



(Al) 



where A^y denote the Laplace operators in the variables x and y, respectively. Because the two 
sides of the equation depend on two different variables, the problem reduces to two ordinary 
differential equations 



and 



A^X(x) 
X(x) 

AyY(y) 
Y(Y) 



= —X-(X + 1), 



(A 2) 



where X • (X + 1) is a separation constant to be determined by implementing the boundary 

(1 2) 

conditions. Equations (A 2), in general, have two linear independent solutions ' \x) and 

fl 2) 

Y\ ' (y). The boundary condition on the tip results in two eigenvalue problems for the separation 
constant X(X + 1): 



and 



y[ 2) (yo) = o. 



(A3) 



The solution of equation (A3) is a countable set of values X^ and X^ with h = l,2, ... and 
k = 1, 2, . . ., which we can order starting from the smaller one. By superposition, the general 



(A 4) 



solution of the Laplace equation fulfiling the boundary condition on the tip reads 

A,, 

+E^ 1) -^ ) ^)+4 2) -xf ) ©]yf ) (y), 

h 

(1 2) 

where the coefficients A K h ' k in general depend on d, &d and are determined by the boundary 
condition on the plane. The conducting plane is generally described by an equation of the type 
x = V(y,d) and does not necessarily coincide with a surface of the type y = yo- By implementing 
the boundary condition on the plane, we obtain the following condition: 



(2) . v(2)/ 



+ EK • x * \ny,d)) + Af ■ xf\V{y,d))]Yf(ii), 



(A 5) 



which must be fulfiled for any y. Consequently, an infinite set of equations for the infinite number 

(1 2) 

of coefficients A h ' k is obtained. 

Scaling hypothesis. The analysis of the boundary condition equation (A 5) determines the scaling 
properties of the potential with d and <P d . We note that the scaling property with & d follows 
directly from the fact that equation (A 5) must hold for any value of d and y, such that each term 
in the sum must scale with 0j as 



(1,2), 



%k 



(A 6) 



The scaling behaviour of the potential with d is determined by the specific geometry of the 
considered problem. 



(a) Conical tip 

Differential equations. We introduce spherical coordinates (r, ft, 0) with respect to the origin of the 
coordinate system [26] (figure 2), such that the surface of the cone is a surface of constant polar 
angle ft = fto. With the separation Ansatz <P{x,y,z) = R(r) • @{ft), we obtain the following Laplace 
equation: 



0 R"(r) R'ir) 



R(r) 
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(A 7) 



the lowest one being the Legendre differential equation. 

General solution. For any value of X, the linearly independent solutions of equation (A 7) are 



R ( l\r) = r\ Rf ) (r)=-1 T 



and 



0 (1) (COS # ) = P X (COS 0 W (COS !? ) = Q X (COS !? ), 



(A 8) 



where Px(cosft) and Qx(cosft) denote the Legendre functions of the first and second kind [27]. 
Some of these functions are, however, not suitable as solutions of our conical tip problem. In fact, 
for non-integer values of X, the Legendre function of the first kind is regular for all angles 0 < 
ft < 7i , but it is singular for ft = tt. However, this singularity is excluded from Q so that Pa (cos ft) 
can be admitted as solution. Legendre functions of the second kind, however, are singular for 
both ft = 0 and ft = tt, such that we must exclude them, as the line ft = 0 is part of Q. Similarly, 



R( 2 )(r) has a singularity at r = 0 and therefore cannot enter the general solution of our problem. 
The general solution of equation (A 7) reads, therefore, 

<P{x,y,z) = J2M®d,d) - r x • P A (cos ft), (A 9) 

where the sum runs over a countable set of values X, which is determined by the boundary 
condition on the cone. 

The eigenvalue problem: boundary condition on the tip. The allowed values for the separation constant 
k are obtained by implementing the boundary condition on the cone 

P x (cosft 0 ) = 0. (A 10) 

The above equation has an infinite countable number of solutions ^(fto), k = 1, 2, 3, . . ., which can 
be ordered with increasing magnitude and have been calculated as a function of fto, for example, 
in ref. [17]. For sharp cones (fto the smallest index X\ can be written approximately as a 
function of the the full angle of aperture of the cone (l>o = 2tt — 2#o 

^ )ss pskw- (A11) 

In summary, the general solution of the electrostatic problem pertaining to a conical tip reads 

<p(r, ft) = ^A k (d, 0 d )r Xk P h (cos ft). (A 12) 

k 

Scaling behaviour: the boundary condition on the plane. The coefficients A k are uniquely determined 
by imposing the boundary condition on the plane placed in front of the tip, defined by the 
equation r cos ft = d with 0 < ft < tt/2: 



COS ft 



P Xk (cosft) = 0 d . (A 13) 



Because equation (A 13) must hold for any value of the parameter d, the ^-dependence must 
cancel out from each term on the left-hand side. This requirement produces the scaling law for 
the parameter d, which combined with equation (A 6) yields 

A k (d,0 d ) = A k (l,l).0 d -d- x K (A 14) 

The potential in the vicinity of the tip apex. The infinite set of constants A k (\, 1) is determined by 
requiring that equation (A 13) is fulfiled for 0 < ft < tt/2. In principle, by multiplying both sides 
by, for example, Px q (cosft) and integrating over the range 0 < ft < tt/2, one can transform this 
boundary condition into a system of an infinite number of coupled linear equations for the 
infinite number of constants A k (\, 1). Unfortunately, finding A k (\, 1) is complicated by the fact that 
the functions of cos ft appearing in the sum do not obey simple orthonormality relations in the 
range 0 < ft < tt/2, so that the diagonalization of the resulting matrix is a cumbersome numerical 
problem which can be solved only approximately. It is, however, not our purpose to find the entire 
set (or even a large number) of the constants A k (l, 1). In fact, of relevance to the discussion of §3 is 
only the potential close to the tip apex and along the tip axis z, which is given by the leading term 
in equation (A 12). Provided that the term equation (A 12) corresponding to X\ is not vanishing, 
we have 

0(z)~0 d -d- kl z kl . (A 15) 

The non-vanishing of this term can be proved using the maximum principle for harmonic 
functions, which states that any harmonic function on a domain Q with boundary d£2 takes its 
maximum and minimum values necessarily on the boundary. This implies that A\ ^ 0: otherwise, 
the leading term would be the one corresponding to Xi> but Px 2 also vanishes again for ft < ft$, 
thus violating the maximum principle for harmonic functions. Note that this expression for the 
potential is valid for any distance d, as, in the vicinity of a cone-singularity, r can be made 
arbitrarily small. 



(b) Sphere on the cone 



17 



This is a slight modification of the conical problem but allows an exact model of a rounding of the 
tip apex, which is expected for realistic tips. 

General solution. The cone acts as a post supporting a small sphere of radius a. The sphere replaces 
the point singularity at the apex of the cone. For the computations in this particular section, it is 
convenient to set the origin of the coordinate system at the centre of the sphere, such that we write 
{x,y,z f = z + a), where z is measured from the tip apex (figure 3). In this coordinate system, the 
position of the plane is given by z p = d + a. The general solution of this axially symmetric problem 
follows the same lines as in the previous subsection. In this case, the general solution reads 

<P(r,&) = J2M^,4>d) (-) -(-) JMcostf), (A16) 

k 

where A^(#o)/ k = 1, 2, 3, . . . are the solutions of equation (A 10). 

(i) Scaling behaviour: distant planes 

The boundary condition on the plane requires 



/ z p \ h _ / a COS ft \ 
\flCOS ft J \ Zp / 



P), k (cos ft) 



(A 17) 



for all 0 < ft < it II. In order to find the scaling behaviour with & d ,z r v , however, it is sufficient to 
consider the boundary condition along the tip axis and thus we set ft = 0 in equation (A 17). For 
z p ^> a, equation (A 17) reduces to 

Y2 A k(z' F »a,<t>d)(j) V(l) = * d , (A 18) 



wherefore the following scaling relation holds 

A k (z f F ^>a,0 d )~0 d . Qjr) • 



(A 19) 



Accordingly, the leading order term of the potential in the vicinity of the tip apex and along the 
tip axis reads 



(A 20) 



When written in terms of the coordinates x,y,z on the tip apex, equation (A 20) recovers 
equation (2.4) exactly, where we have exploited the fact that z p ~ d, which holds for distant planes. 

(ii) Scaling behaviour: near planes. 

As described above, it is sufficient to set ft = 0 in equation (A 17). For planes close to the tip, 
(z p — a) /a <^ 1, wherefore 



C P ) (4) 



to leading order. From equation (A 21), we deduce the scaling behaviour 
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(A 22) 



The leading order term of the potential along the z ; -axis and in the vicinity of the tip apex is thus 
given by 

*(zf)~* d .^^, (A 23) 



b 
(3 



which is exactly equation (2.4) in terms of the coordinates (x,y,z). 



(c) Hyperboloid of revolution 

Differential equations and general solution. The tip surface is modelled as a hyperboloid of 
revolution, obtained by rotating a hyperbola around the axis containing its focus with coordinates 
x = y = 0 and z' — —a (figure 3). In order to implement the boundary condition on the tip surface, 
it is convenient to introduce the prolate-spheroidal coordinates (/x, y,0) [28], with /x e [0, oo[, 
v e [0, 7r], and 0 e [0, 7jz\, defined through 



and 



x = a sinh \x sin v cos 0, 
y = a sinh /x sin y sin 0 
z' = 0 cosh \i cosy. 



(A 24) 



Surfaces of constant \x are ellipsoids of revolution with foci at z' = ±a, where the limiting case 
\x = 0 is a straight line connecting the foci. Surfaces of constant v are hyperboloids with focus 
at +a for v e [0,7i/2] and at —a for v e [tt/2,7t]. In particular, the surface of the tip is v = vq. The 
corresponding Laplace equation reads 



a 2 ^ 
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dr, (§2 _ i)(i _ ,2) , 



$ =0, 



(A 25) 



where we have introduced £ = cosh/x and = cosy. The rotational symmetry of the problem 
implies that the potential is independent of the angle 0. With the separation Ansatz @(%,r]) = 
S (i-)N(rj), we obtain the following equations: 



and 

with solution 



(i-r)g^ sG)]+x(x + i)SG) = o 



A 



(A 26) 



(A 27) 



where the sum extends over suitable values of X, determined by implementing the boundary 
condition on the tip. 

(i) Scaling behaviour: distant planes (z f ? /a » 1) 

In the limit of distant planes, we must exclude terms proportional to Qa(£) in equation (A 27), as 
they diverge along the line y = 0 (x = y = 0,a < z f < oo). The boundary condition P^(cos vq) = 0 
is satisfied by the same values that we have obtained when considering conical tips. The 
boundary condition on a distant plane involves finding the behaviour of Px k and Qx k for large 
arguments. For this purpose, we let \x -> oo and a -> 0 by keeping r = a • cosh \x constant. To 
visualize the significance of this limit, consider that in this limit the foci are 'moved' towards the 
origin and a cosh \x « a sinh fi, so that constant-^ ellipsoids become constant-r spheres. Formally, 
in the limit \i -> oo with r = a%= constant, prolate spheroidal coordinates are in fact the same 
as spherical coordinates. To appreciate this point, we rewrite the differential equation for £ in 
terms of r 



4(((D 2 - 1 )s a - w )- ita > w=a 



In the limit a — >► 0, equation (A 28) simplifies to 

r*El(r) + 2rE' Xt {r)-X k E kk (r) = Q, 

whose solution reads 



(A 28) 

(A 29) 
(A 30) 



For distant planes, S^ 2 \zp/ a cos v) becomes negligible and the boundary condition on the plane 
simplifies to 

*- = E^ff(^)-^(«»v), (A31) 

which, in virtue of the asymptotic behaviour of 3^(r), requires the scaling property 

A Xk (0 d ,z f F )~0 d .^ \ (A 32) 

Potential close to the apex. The potential in the vicinity of the tip apex does not contain Q^(cosh /x) 
which is singular in the interval z' e [—a, a], corresponding to \x = 0. In addition, we point out 
that in the vicinity of cos vo both P^(cos y) and Q^(cos v) are analytic functions having a Taylor 
expansion starting with a linear term in cos v — cos vq, z f 0 =a cos vo being the z'-coordinate of the 
tip apex. Thus, to leading order in {z! — z^/a, the potential in the vicinity of the hyperboloid apex 
and along the hyperboloid axis reads 

* w ~*'-un^)- (a33) 

This is equivalent to equation (2.4) in terms of the x, y, z-coordinates (recall that z f = z + a cos vq). 



(ii) Scaling behaviour: near plan [\z f ? — z f Q \/a <^ 1) 

In this case, we must exclude terms proportional to Qa(m) i n equation (A 27), as they diverge 
along the axis \x = 0. The boundary condition along the surface of the hyperboloid selects indices 
kh and for which either P^(cos vq) = 0 (and Bx h = 0) or Q^(cos vq) = 0 (and A\ k = 0) holds. The 
general solution then reads 

*($,!/) = Y J Ax h PxM ■ Pa*0?) + ^BfcPfcfc) • Qx k M- (A 34) 

We consider the boundary condition on the plane in the case \x = 0 and obtain 

^ = J2 A MV • Ph ( Z i) + ^B^(l) • Qx h (A 35) 

Because for Zp close to the tip apex both Px h and Qx k can be approximated by linear terms in 
(z p — Zq)/0, equation (A 35) implies the following scaling relations: 

and B ^^7T7- ( A36 ) 

Z p Zq Zp Zq 

Potential near the apex. To leading order, the potential in the vicinity of the hyperboloid apex and 
along the hyperboloid axis reads 

0(2^^.^ (A 37) 

Zp Zq 

which is equivalent to equation (2.4) in terms of the coordinates x,y and z. 



(iii) Scaling behaviour: confocal plane (Zp = 0) 

This case is special, because the conducting plane coincides with the hyperboloid defined by 
v = tt/2, and thus it is confocal to the tip surface. Boundary conditions are therefore homogeneous 
in the variable v, such that the solution does not depend on jjl. This requirement determines 



uniquely that A = 0 and Bq = 0 is the only possible solution. Accordingly, the general solution 
simplifies to 

C 0 P 0 fa) + DqQoW = C 0 + ^ In f ! + C ° SV \ (A 38) 

The constants Co and Dq must be chosen in order to satisfy the boundary conditions <P{v = vq) = 0 
and 0(v = tt/2) = +0d, so that the potential reads, exactly [29] 



0(/i,cosv,(P) = 0 d - — — ^— -lnf „ +cosv Y (A39) 

\ 1 — cos V ) 



ln((l + cos v 0 )/(l - cos vq)) 



Close to Zq = a cos vq and along the z'-axis we have, to leading order in z' — z r Q 



ln((l + cos vq)/(1 — cos vq)) a • (1 — cos 2 vq) 



^•i^kr- (A40) 

Here, we have used the relation a 2 cos 2 vo = d 2 and cos vo = y/d/(d + Ro)/ Ro^d being the radius 
of curvature of the hyperboloid. This result shows that the linear dependence on 1/d of the 
potential crosses over to a logarithmic dependence when the planar electrode crosses the confocal 
plane. This is consistent with the idea that equation (1.1) corresponds to a logarithmic dependence 
on d when X = 0, which we put forward when discussing the paraboloid solution for d^>Ro. 
When the plane becomes a distant one a true power law intervenes. 

A remark. One might be tempted to use the method just illustrated to solve the electrostatic 
problem of an infinitely thin semi-infinite needle placed vertical to a conducting plane [30]. In fact, 
a semi-infinite needle can be viewed as the limit of an axially symmetric hyperboloidal tip with 
angle of aperture coo ~> 0- We point out that in the hyperboloidal problem, X\ becomes arbitrarily 
small when coq -> 0, both being well-defined and finite as long as coq ^ 0. However, for coo = 0/ the 
boundary condition on the tip cannot be fulfiled, as the equation P^costt) = 0 has no solution 
(in particular, k\ = 0 is not a solution of this equation). This result shows that it is impossible to 
run a junction between a truly semi-infinite one-dimensional wire and a planar counterelectrode at 
a finite potential difference, at least in a three-dimensional space. 



b 



CD 
CD 



(d) Paraboloid of revolution 

Differential equations and general solution. We now model the microtip as paraboloid of revolution 
by introducing three-dimensional parabolic coordinates [31] (figure 4) 



and 



X — UV COS0, 

y = uv sin0 

2 2 

, vr — ir 
z — — - — / 



(A 41) 



where u, v e [0, oo[ and cp e [0,27r[. Surfaces of constant u are paraboloids of revolution with 
negative curvature, the limiting case u = 0 corresponding to the negative z-axis. Surfaces of 
constant v are paraboloids of revolution with positive curvature, where the limiting case v = 0 
corresponds to the negative z-axis. The microtip is taken to be a paraboloidal surface of constant 
u = Uo or, equivalently, it is described by the equation z' = — (l/2)(r/Ro) + ^o/2, where Ro = Uq 
is the radius of curvature of the paraboloid. The coordinates of the tip apex are x = y = 0 and 
= uq /2. The Laplace equation in these coordinates reads 



1 



l d 

U du 



du 



-H't) 

v av \ dv J 



,0=0. 



(A 42) 



By exploiting the rotationally symmetric Ansatz @(u, v,(p) = U(u) • V(v), we obtain the following 
differential equations: 

u 2 U"(u) + uU\u) - n 2 u 2 U(u) = 0 
v 2 V"(v) + vV'(v) + n 2 v 2 V(v) = 0, 



and 



(A 43) 



where n 2 is the separation constant. The solution of equation (A 43) can be written in terms of 
Bessel functions Jo(nu) and Yq(hu) and modified Bessel functions Io(nv) and Kq(hv)) of order 0 



0(u, v) = Yy A nM n ku)h{n k v) + B mk Y 0 (m k u)I 0 (m k v)] 

k 

and 0 d = A n k h (nu^z'-p + v 2 ^ h(n k v) + B m Y 0 (m k ^2z r F + v 2 ^j 



I 0 (m k v) 



(A 44) 



where n k and m k are defined by Jo(n k uo) = 0 and Yo(m k uo) = 0, as a consequence of implementing 
the boundary condition @(uo, v) = 0 on the tip surface. To obtain equation (A 44), we have 
implemented the boundary condition on the plane (u 2 — v 2 )/2 = Zp as well. 

(i) Scaling behaviour: near planes (y^Zp — uq > uq). 

We consider the boundary condition along the line v = 0 and exploit the fact that for near planes, 
^2z p is close to uq. Because the functions Jo(n k u) and Yo(m k u) vanish at Uq, we have to lowest 
order 



Jo (n k ^2^j - 
(m k ^2z^j 



and Yq 
These equations imply the scaling behaviour 



/2z p - Uq 



(A 45) 



An k ,B nk ' 



'*d-r 



(A 46) 



and the potential close to the tip apex and along the z^axis reads, therefore, 



(A 47) 



In terms of the x,y and z-coordinates, equation (A 47) recovers equation (2.5) exactly. 

(ii) Scaling behaviour: distant planes (yftep — uq » uq) 

The formulation of a scaling law in the case of a distant planar counterelectrode is difficult, 
because Jo and Yq share the same asymptotic behaviour and are both oscillating, so that 
cancellations might occur. However, a strategy which foresees to approximate a large but finite 
distant plane with the flat portion of a confocal paraboloid u = u\ = ^2z p might give a useful 
approximation. This specific choice of the counterelectrode geometry has the advantage that 
boundary conditions are homogeneous in the coordinate u, such that the solution solely depends 
on the variable u. Under these circumstances, equation (A 42) simplifies to 



u z <P"(u) + u<P'(u) = 0. 



(A 48) 



By accounting for the boundary conditions on the tip and on the parabolic counterelectrode, the 
solution of equation (A 48) reads 



*(«) 



ln(Mi/w 0 ) 



•In 



(-)■ 



(A 49) 



Close to uq and along the tip axis, the potential reads, therefore, 



•tf-'-isU)-^ (A50) 



which is equivalent to equation (2.5) in terms of the coordinates x,y and z. 



(e) Cuspidal tip 

Differential equations. In this case, it is convenient to introduce bi-spherical coordinates {a, p, 0) 
(figure 2), defined through [32] 

a sin a cos 0 

x = 



cosh p — cos a 

a sin a sin 0 

^ cosh /3 — cos a 

. a sinh p 

z = , 

cosh /3 — cos a 

where 0 <a <tt, — oo</3<oo, 0<(p <2tt and a is the focal distance. Surfaces of constant a are 
tori with radius a/ sin a and centre on the z! = 0-plane at a distance a cot a from the origin of the 
coordinate system 

(J x2 +V 1 - a cot a) + z /2 =— (A 51) 
V v / sin a 

whereas surfaces of constant /3 are spheres of radius a/ sinh /3 centred at (0, 0, a coth /3) 

a 2 

(t! - a coth p) 1 + {x 2 +y 2 )= — . (A 52) 

sinh p 

In these coordinates, the tip surface is given by a = ao > 0, such that the associated boundary 
condition is <i>(ao,/3,0) = 0. The corresponding Laplace equation is separable if the following 
rotationally symmetric Ansatz is used 



<P(a, P) = y^cosh^ -2 cos a • A(a) • B(fi). 
The functions A(a) and B(P) satisfy the following equations: 

1 d / 1 cL4\ t ^ A 

+ k(k + 1)A = 0 



sin a da \sina da 

d 2 B / 1\ 2 „ 
and = ( A + 2 ) 



d^ 



(A 53) 



with solution 



&(a,P) = y/2coshp -2 cos a 

x J][A A P A (-cosa) + B A Q A (cosa)] 

A, 

x [C X e^ +1 ^ + D k e-^ +l ^l (A 54) 

Divergences at a = tt {—a<z'<a) are avoided by setting Bx = 0. 

Boundary condition on the tip. By implementing the boundary condition on the tip, we obtain the 

equation Px(— cosao) = 0, with solutions These indices coincide with those computed 

for the conical tip model with aperture angle coo — 2ao- 

Scaling properties. Let us consider a plane at z! — Zp close to tip apex. Note that this plane cannot 
extend indefinitely in the xy-direction, because it would encounter the boundary represented by 
the tip itself, which is kept at a different potential. Relevant in this case is the behaviour of the 



solution for ft -> — oo. Under these circumstances, sinh(/3) ~ e & /2 and coth(/3) ~ — 1, such that 
equation (A 52) simplifies to 

(x 2 +y 2 ) + (z^+fl) 2 = p^ / (A 55) 

which is the equation of a sphere with radius 2a e@ centred at the tip apex. Note that for the above 
approximation to hold, z p needs to lie in the interval [—0,0]. Along the tip axis, equation (A 55) 
implies /3 « ln((z p + a) /2a), such that we must set D\ = 0 in equation (A 54) to avoid divergence 
at z p -> —a. The boundary condition along the axis reads accordingly 



00 / z f +a\ Xk 

®d = Yl Cx k (<P d/ z f F )P Xk (l) , (A 56) 

k=i \ a J 



wherefore 



^—'dzTiY- (A57) 

Finally, the leading order term of the potential along the tip axis and in the vicinity of the apex 
reads 

which is equivalent to equation (2.2) in terms of x,y and z (recall: z' p + a = A, figure 2). 

Appendix B. The tunnelling current density 

In this appendix, we present the details regarding the computation of the integrals appearing in 
G(E f ) andJo(Ef). 



(a) Non-analytic potential 

The potential within the tunnelling barrier writes, to leading order, <t> (z) = Ax 1 • @d ' (z/d) Xl , where 
Ax 1 is a dimensionless constant. Referring to figure 5, we set z\ — 0 and 



Z2=d 



A x 



By introducing the dimensionless variable 

y= A—) ' 



we obtain 



G(E f ) = 2 



87t 2 m 



l87T 2 m(p 



lAi 



lAi 



-r(i + iAi) 
2 - r(3/2 + 1/A.i) 



:27T 3 / 2 



r(i + iAi) 
r(3/2 + iAi) 



A ( 



-lAi 
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(B3) 



and, similarly, 



7o(Ef) = 



47vem I 87T 2 m d • ^lAi-l/2 

47rem ( j87T 2 m d'(p 1/Xl ~ 1/2 
h 3 



h 2 (A kl e* d )V* r(l/2 + lAi) 




• rg/2 + iAi) 

.v^ni+i/Xi) 



2/A.i 



(B4) 



We observe that both the expressions for G(£f) and /rj(Ef) depend on the scaling variable 
d~ Xl , whose natural unit turns out to be <p A^ Xl /e. 

(b) Rounded tips: linear potential 

In this case, the potential is given by <P{z) = Ax^a/d)^ 1 ■ 0 d - z/a for distant planes providing the 
boundary condition and by &(z) = A\ • 0 d • z/d for near planes. Therefore, we set z\ = 0 and 

d - (p 



Z2 : 



and 



Z2 = 



a - cp 



(near plane) 

(distant plane). 



(B5) 



A Xl e<P d (a/d)^ 

The computation of G(£f) and /o(Ef) reduces to the computation of the elementary integrals 

rl j 



1 — y 1 I dy and 



dy, 



(B6) 



such that one obtains 



and 



G near (0 d ,d) = 

G distant (<P d ,d) = 
J^(0 d/ d) = 



8tt 



A 1 e0 d • Ayld 

e • cp (A 1 e0 d • Ay/d) 2 



8jt h-A 2 



8tt a 
T " A 
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cp A Xl e0 d • (a/d) x 



e-(p A% (A Xl e0 d -(a/d)^) 2 
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(c) The appearance of A v 

A most remarkable feature common to all equations in this appendix is the appearance of a 
characteristic length A<p, which is the De Broglie wavelength associated with an electron with 
an energy corresponding to the potential energy barrier (p. The simplest way to illustrate this 
somewhat surprising appearance of A v is to consider the Gamov factor for the elementary 
triangular barrier in the planar geometry [33] 

V2m (p 3/2 
3 e-h ~T' 

where F is the electric field and the superscript el used in ref. [33] probably stands for 'elementary 7 . 
In the planar geometry, one has F = 0 d /d. Accordingly, G el writes 



(B8) 



b 
(3 



Note that the Gamov factor must be dimensionless, as it is an exponent. To rewrite G el as a 
manifestly dimensionless quantity proceed as follows: 



(i) (p 3 / 2 is written as cp • cp 1 / 2 . 

(ii) The first factor (p has the units of energy and is used to cancel the dimension of e<S>&- 

(iii) The second factor (p 1 / 2 is associated with the remaining terms to create a quantity with 
the dimensions of an inverse length, which cancels the dimension of d. One finds that this 
quantity is exactly the inverse of the characteristic length A<p. 

This shows that A<p was just hidden in the widely known expression for G el . Having (proudly) 
demonstrated the appearance of A v in the Gamov exponent, we discovered that the same length 
was already introduced into a Gamov exponent more than 80 years ago by J. R. Oppenheimer 
(eqn (1) in ref. [12])! 
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